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We investigate the influence of columnar defects in layered superconductors on the thermally 
activated penetration of pancake vortices through the surface barrier. Columnar defects, located 
near the surface, facilitate penetration of vortices through the surface barrier, by creating "weak 
spots", through which pancakes can penetrate into the superconductor. Penetration of a pancake 
mediated by an isolated column, located near the surface, is a two-stage process involving hopping 
from the surface to the column and the detachment from the column into the bulk; each stage is 
controlled by its own activation barrier. The resulting effective energy is equal to the maximum of 
those two barriers. For a given external field there exists an optimum location of the column for 
which the barriers for the both processes are equal and the reduction of the effective penetration 
barrier is maximal. At high fields the effective penetration field is approximately two times smaller 
than in unirradiated samples. We also estimate the suppression of the effective penetration field 
by column clusters. This mechanism provides further reduction of the penetration field at low 
temperatures. 



I. INTRODUCTION 

The properties of the Abrikosov vortex state of type II 
superconductors with artificially manufactured columnar 
defects attract a great deal of current attention. Moti- 
vated originally by a technological quest for enhanced 
vortex pinning, superconductors with columnar defects 
revealed a vast diversity of remarkable phenomena. The 
possibility to introduce controlled disorder and to tune 
parameters (such as vortex density and interactions be- 
tween vortices) has made them one of the favorite exper- 
imental systems for studies of the statistical mechanics 
and dynamics of a glassy state (see recent review article 
|l| and references therein) . 

On the other hand, the magnetic response of high- 
temperature superconductors is known to be controlled 
to a large extent by the creep, of vortices over the Bean- 
Livingston surface barrier,0cl an important manifesta- 
tion of which is the exponential temperature dependence 
of the effective penetration field. It was recently observed 
that columnar defects can strongly facilitate the creep of 
pancake vortices ever the surface barrier and reduce the 
penetration field.Q Usually pinning by disorder and sur- 
face barrier are considered to be competing effects that 
alternatively control the magnetic response in the vor- 
tex state. In this paper we consider an interplay be- 
tween surface and bulk pinning and develop a theory for 
the disorder-assisted surface creep in highly anisotropic 
superconductors focusing on the case of randomly dis- 
tributed columnar defects. This effect is somewhat anal- 
ogous to tunnelling of auantum particle in the presence 
of sub-barrier disorder B 

It is clear from a general consideration that surface im- 
perfections create weak spots, facilitating penetration of 
vortices through the Bean-Livingston barrier. However, 
a quantitative theory for imperfection of arbitrary kind is 



not available. Special kinds of surface irregularities have 
been considered in Refs. p|-pT|. We address well defined 
surface disorder created by controlled irradiation with 
heavy ions. In this case vortices enter the sample near 
the weak spots where columnar defects are located close 
enough to the surface to suppress the surface barrier. 
Vortex penetration consists of two steps (see Fig. [I]): (i) 
the capturing of the vortex onto a near-surface column 
and (ii) the detachment of the trapped vortex into the 
bulk . The resulting effective barrier for pancake pene- 
tration via an isolated column is the maximal value of the 
two barriers corresponding to the above processes. The 
capturing process has, in its turn, a two-channel char- 
acter and may occur either via the direct motion of a 
pancake to a column or via the nucleation of an antivor- 
tex at the column and its subsequent advance towards the 
surface (see Fig. Q). The capturing process is controlled 
by the channel with the smallest barrier. If the mediating 
column is located far from the surface, the penetration 
process is controlled by the transfer of a vortex from the 
surface to this column. When the column is sufficiently 
close to the surface, the controlling barrier corresponds 
to detaching a vortex from the column to the bulk. For 
every external field there exists an optimal location of 
the column for which the reduction of the barrier is max- 
imal, and penetration of pancakes into the sample occurs 
mainly via such optimally placed columns. 

One can expect that the surface barrier may be sup- 
pressed even more efficiently in very rare places where 
several columns happen to be near the surface. The net 
contribution to the penetration rate from such events is 
determined as a balance between their small probabil- 
ity and the strong loca l suppression of the barrier. We 
estimate in section III the collective suppression by the 



column clusters. At low temperatures the collective sup- 
pression always becomes more efficient than the single- 
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column mechanism. 
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FIG. 1. Mechanisms of penetration of pancake vortex from 
the surface into the bulk of a superconductor: (a) Nucleation 
of the vortex at the column via the motion of vortex from the 
surface (b)Nucleation of vortex at the column via motion of 
an antivortex from the column to the surface (c) motion of 
the vortex from the column into the bulk 



II. VORTEX ENERGY NEAR SURFACE IN 
PRESENCE OF AN ISOLATED COLUMNAR 
DEFECT 

We consider an irradiated superconductor with insulat- 
ing columnar defects oriented along the c-axis in external 
magnetic field also applied along the c-axis. Let an iso- 
lated columnar defect have radius R and be located at a 
distance L from the surface of superconductor (see Fig. 
||) . The energy of a pancake vortex located between the 
surface and column at distance x from the surface con- 
sists of two parts: the direct interactions with the col- 
umn and the surface, Ui nt {x), and the interaction with 
the Meissner current, Uj(x). We introduce dimensionless 
variables measuring length in units of i?, energy in units 
of sen = s$q/ (47rA) 2 , current in units of R , and 

magnetic field in units of . $ ° . Here A is the in-plane 
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London penetration depth and s is the interlayer spacing. 
We also use the notation I = L/R. 



FIG. 2. Geometry: a pancake vortex between the surface 
and a columnar defect. Interaction with the surface and the 
column can be described in terms of an infinite set of positive 
and negative images obtained by multiple subsequent reflec- 
tions by the surface and the column (only the first two images 
inside the column are shown). 



A. Direct interaction with surface and column 

The interaction of the vortex with either the plain sur- 
face or with an isolated column can be calculated by 
virtue of the image technique!! 2 ! The interaction with the 
surface is obtained by placing a negative vortex at the 
point — x and the interaction with the isolated column 
can be obtained by placing a negative vortex at point 
I — -r^— inside the column and a positive vortex at the 

n 

column center. Lj In the case where the vortex is con- 
fined between the column cavity and the surface, adding 
these images only would not solve the problem because 
the currents due to the surface image do not satisfy the 
boundary conditions at the surface of the columnar cav- 
ity and vice versa. To compensate the currents of the 
surface image, one has to add the its image (reflection) 
inside the column. This eliminates vortex in the column 
center and adds a vortex at the point / — . After that 
we have to add the surface image of this vortex at the 
point —I + jqL. Continuing these reflections, we obtain 
an infinite set of positive and negative images. We label 
coordinates of the positive (negative) images inside the 
column after the n-th double reflection as u+ (u~) (every 
such image has a corresponding surface image of opposite 
sign at -u^(-ti^)). The n + 1-st image is obtained by 
reflecting the n-th image with respect to the surface, and 
then reflecting it again with respect to the column. This 
yields the following recurrency relations 

u n+l = I ~ 1~. - j U n+1 = ' — 7~ + ) (1) 

with 

u I = 1 - 7^—; uf = 1- — — , 

I — x I + x 

As follows from ([j]), both the negative and positive im- 
ages approach the same limiting position = — 1 



2 



as n — > oo. The interaction energy is then expressed as 
an infinite series in the coordinates of all images: 



U int (x) = In 
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(u n + x) (u£ - x) 



Note that this expression is valid for a vortex located at 
either side of the column, i.e., for both < x < I — 1 
and x > I + 1. The supercurrent distribution around 
the column containing the trapped vortex coincides with 
the distribution corresponding to a vortex placed at the 
point \Jl 2 — 1 and its surface image. The energy of this 
state is given by 



£7^(0 = ln(l + V* 2 - 1) 

The energy, corresponding to an antivortex at the point 
x, which controls the competing process of antivortex 
motion from the column to the surface (we remind that 
in this process a vortex-antivortex pair nucleates at the 
column, then the antivortex leaves the column and hops 
to the surface, while the vortex remains trapped in the 
column) is obtained analogously and reads 



U%\x) = U int (x) + 2ln ^=^ 



+ ln(J + Vl 2 ~ !)• 
(3) 



Here the second term in the RHS describes the interac- 
tion of antivortex with the trapped vortex while the third 
term gives the self-energy of the trapped vortex. 

To solve the recurrency relations (yj), we introduce a 
new variable /„ for both the positive and negative images 



u n = ut as 



fn 



Vi 2 - i - u„ 

and transform Eq. (jl|) into a simple relation 

fn+l — 0? fn 



with a = I + \/l 2 — 1 > 1, which can be easily solved: 
/„ = a 2n fi. This allows us to obtain closed analytical 
formulas for the image coordinates: 



= sinh b 



= sinh b 



x + tanh(nfe) sinh b 
tanh(n&)x + sinh b 
x — tanh(nfo) sinh b 
tanh(n6)x — sinh b 



(4) 



with b = In (l + \/l 2 — l) and I = cosh b. Taken together, 
Eqs.(||) and (|J) solve the problem of finding the energy 
of a vortex located on the perpendicular to the surface 
passing through the center of the column. This result 



can also be obtained using complex plane representation 
and the conformal mapping 



= In 
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It transforms the semispace x > with a circu- 
lar hole at z = (1,0) into the rectangular area 
{0 < wi < b, — 7r < W2 < 7r} with the periodical bound- 
ary condition along the u>2 direction. 



B. Interaction with the Meissner currents 

A column (cylindrical cavity) placed near the surface 
disturbs the pattern of the screening supercurrent in- 
duced by the external magnetic field and changes accord- 
ingly the contribution to the vortex energy arising from 
surface screening current. The current j(r) has to satisfy 
diyj = 0. This means that it can be expressed in terms 
of a supercurrent potential, ^j(r), as 



dx 



3x 



dy 



We consider the situation where all relevant distances are 
smaller than the London penetration depth so that we 
can neglect screening effects and assume curlj = 0, which 
implies that the potential satisfies the Laplace equation 
A(f>j — 0. The problem of finding the current distribu- 
tion is thus equivalent to the problem of the polarization 
of a metallic cylinder near a metallic surface-,by the ex- 
ternal electric field parallel to the surface.EJ Since the 
normal components of j should vanish at the surface and 
at the column, both the external boundary of the su- 
perconductor and the boundary of the column should 
be equipotential surfaces. Setting (pj(0,y) = 0, we de- 
fine (f>j(r) as the interaction energy of a vortex at the 
point r with the Meissner current. Note that at large dis- 
tances from the column, the x-component of the current 
should vanish, thus, far from the column j = (0, j) with 
j = cHj (4ttA). The current distribution near an isolated 
cylinder can be obtained by putting vortex dipole (j, 0) 
at the center of the cylinder, this dipole induces poten- 



tial <f>j — —jx (l 



To satisfy the boundary 



condition at the surface one has to add surface image of 
this dipole, i.e. to put dipole (—j, 0) at x = —L. How- 
ever, the currents of this surface image do not satisfy the 
boundary conditions for the column and we again have 
to add the column image of the surface image. Continu- 
ing this process, we again obtain an infinite set of dipole 
images inside the column. Note that this sequential re- 
flection in the column does not preserve the magnitude of 
a dipole. The reflecting pair of opposite vortices located 
near the point (x, 0), we derive that the magnitude of col- 
umn reflection of dipole is smaller by factor 1/ (I — x) 
than the magnitude of the original dipole. Denoting the 
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coordinate of the dipole resulting from the n double re- 
flections and its magnitude as jp n , we derive the 
recurrency relations 

1 ,1 

(l + x n y l + Xn 



The vortex-Meissner current interaction energy be- 
comes (f>j(l — 1) = —jy/P — 1, when the vortex is trapped 
by the column, and the total energy of the trapped vortex 
reads: 



U tr (l,j) = ln(l + y/p - 1) - jVl 2 ~ 1- 



with x\ = I, pi — 1. Analytical solutions to these equa- 
tions are given by 



x n = \]l 2 — 1 cothnfr, 
I 2 -I 

Pn 



sinh nb 



(6) 



where, again, b = In (Z + yW—1) . The potential can be 
represented as an infinite series 



n=l 



[X Xn) p n 



(x + x n )p n 



(x-x n f+y 2 {x + x n f + y 2 



The interaction energy of the vortex, located at the line 
y = 0, with the Meissner current is given by 



Pn 



<t>3 ( X ) = §3 (»> 0) = -jx + 2xj ^2 



Fig. H illustrates the current distribution for 1 = 2. 
When column is located sufficiently close to the surface, 
I — * 1, the current at the surface diverges as ^(0,0) = 
V2j/\/l ~ 1, while the current at the opposite side of 
the column approaches a universal value, j y (l + 1,0) — > 

j(l + 2Er =1 #^)=(- 2 /4)^2.4 7j . 
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FIG. 3. Example of the current distribution for I = 2. 



C. Activation barriers 



As we have already mentioned, the penetration of a 
pancake into an irradiated superconductor happens in 
two steps: hopping from the surface to a column, and de- 
tachment from the column into the bulk. The resulting 
effective barrier is equal to the maximum of the two bar- 
riers for the two processes. In addition, the nucleation of 
a vortex at the column (the first step) can occur via two 
channels: (i) as motion of a vortex from the surface to a 
column or (ii) by nucleation at and the subsequent mo- 
tion of an antivortex from the column to the surface. The 
effective barrier for this first step is, thus, the smaller one: 
the process goes through the easier channel. For every 
fixed position of the column there exists a certain value, 
jcd(l), of the surface current above which the state with 
the one flux quantum trapped in the column becomes 
energetically favorable: 



3cd{l) = 



ln(Z + yp^i) 



This current is plotted in Fig. ^. The magnetic field cor- 
responding to this current is given by 
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FIG. 4. Lines in the j-l plane at which the penetration 
mechanism changes: at currents j > j c d(l) it is energetically 
favorable to put one flux quantum on the column, at cur- 
rents j > jav(l) nucleation of the flux quantum at the column 
occurs by motion of an antivortex from the column to the sur- 
face, the line l op t(j) gives the optimum position of the column 
corresponding to the maximum suppression of the barrier. 
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H cd (L) = 



4vrAVi 2 - R 2 
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L + VL 2 - R 2 
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Using the total energy of the vortex located at the 
point (a;, 0), 

U v (x) = U int (x) + (pj(x), 

one can calculate the activation barrier for the vortex to 
hop from a surface to the column (process 1 in Fig. |l|), 

UAl,j,R/ti)= m ax U v (x). 

Analogously, making use of the expression for the energy 
of the couple, the antivortex at the point (x, 0) and the 
vortex inside the column, 

u av (x) = u£:\ x ) - Mx) - jVp~T, 

one obtains the barrier for the motion of the antivortex 
from the column to the surface (process 2 in Fig. [I]) , 

U2(hj,R/Q= max U av (x). 

0<x<l — l 

The trajectory of the process, i.e. the channel that will 
actually govern vortex penetration, depends on the mag- 
nitude of the applied current. At small currents the 
surface-to-column process dominates while at sufficiently 
large currents the vortex-antivortex mechanism comes 
into play. The characteristic current, 3 av (l), separating 
these two regimes depends on the position of the col- 
umn and is determined by the solution of the equation 
U\(l,j a v) = U 2 (l,3av)- The plot of j av (l) is shown in Fig. 
H together with j c d(l)- The barrier to activate the vortex 
from the column into the bulk of superconductor (process 
3 in Fig. [I]) is given by 



U 3 (l,j,R/0-- 



max x>i+ i \U v (x)\ , j < j c d(l) 
max 2 . >/+1 [U v (x)]-U tr (l,j),j>j cd (l) 



Finally, the total barrier corresponding the channel 
surfaces-column— >bulk reads: 



j, R/0 = max [min («i,« 2 ) Mz] 



(7) 



In the following we will calculate the effective reduction 
of the surface barrier by the column 

5U(l,j)=U(l,j,R/0-UoU,R/0, (8) 
which does not depend on the ratio R/£. Here 
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U {j,R/0=\n- 
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(9) 
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FIG. 5. The evolution of the energy profiles with increas- 
ing screening surface current j for 1 — 2. On the lefthand 
side: solid lines represent the energy profiles, U v (x), for the 
vortex moving from the surface to the column, and dashed 
lines represent the energy profiles for an antivortex, U av (x), 
moving in the opposite direction. 
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FIG. 6. Dependence of the energy barriers on relative dis- 
tance from the surface I = L/ R for the three processes shown 
in Fig. hi for several values of j. 
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FIG. 7. 1-dependence of the total barrier for penetration 
from the surface into the bulk for several values of j. 
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FIG. 8. j- dependence of the energy barrier reduction for 
penetration of a vortex from the surface into the bulk for the 
optimum column location, i.e. for the minimum positions of 
5U(l,j) plotted in Fig. (?]. One can expect that this plot de- 
termines the reduction of the surface barrier in real irradiated 
samples in the single column regime. 

We will now turn to a numerical evaluation of the 
barriers. Fig. ^| illustrates the evolution of the energy 
profiles U v (x) and U av (x) with increasing current for 
1 = 2. As one can see, at small current the position 
of the maximum energy is located at the righthand side 
of the column while at large current it is located be- 
tween the surface and column. Fig. |] shows examples of 
the / dependence of the column-induced barrier changes 
8U\(l,j), SUi{l,j) and SU 3 (l,j) for different j. The bar- 
riers Ui(l,j) and U2(l, j) are rapidly suppressed when 
the column approaches the surface, mainly because of 
the divergence of the surface current. On the other 
hand, the barrier 5lA^(l,j) increases slowly with the de- 
crease of /. Thus, for large / the total barrier is deter- 
mined by 1A\(l,j) (or 142(1, j)) and for small I it is de- 
termined by Uz(l, j). Fig. shows the /-dependence of 
the total barrier 14(1, j) for different I. For every value 
of the surface current (external field) there is an opti- 
mum location of the column l opt for which suppression 
of the barrier is maximum. It corresponds to a tran- 
sition between the two mechanisms of penetration, i.e., 
Uz(l op uj) = miii(Ui(l pt,j),U 2 (loptJ))- The line l op t(j) 



is shown at Fig. || together with j c d(l) and j av (l). Fig. || 
shows the current dependence of the barrier suppression 
for an optimally located column 5U opt (j) = Sl4(l op t,j)- 
This plot represents the main result of the paper. The 
important feature of this dependence is that at high cur- 
rents it saturates at « —0.76. 

An important manifestation of the surface barrier in 
high-T c superconductorSjis the enhancement of the ef- 
fective penetration fieldB When the penetration field 
is limited by thermal penetration through the surface 
barrier, the effective penetration field is determined by 
the equation Uo(H p o) = CT, where the numerical con- 
stant C ~ 20 — 40 is determined be the.jexperimental 
time scale and by the attempt frequency.Btj This gives 
the exponential temperature dependence of the penetra- 
tion field. In irradiated samples the penetration field is 
expected to be shifted to lower value, which is deter- 
mined by Uo(Hp) + SU(H p ) = CT. In the field range 
H p ^> where the barrier reduction approaches the 

constant value —0.76 se<j we arrive at the very simple 
result Hp ps 0A7H p0 . 



III. COLLECTIVE SUPPRESSION OF THE 
SURFACE BARRIER BY COLUMN CLUSTERS 



In the previous sections we considered suppression of 
the surface barrier due to an isolated columnar defect 
near the surface. Such a mechanism does not give a 
full picture. The barrier may be suppressed even more 
substantially in small number of places where several 
columns occur to be close to the surface. The resulting 
enhancement of the average penetration rate comes as a 
trade between the suppression of the barrier and small 
probability of such event. In this section we estimate the 
contribution to the penetration rate and the reduction of 
the effective penetration field caused by clusters of de- 
fects. The optimal gate for the vortex entry will appear 
as several columns line up next to each other forming 
thus a cut thrusting the sample normal to the surface. 
The width of the cut is 2R and its length is 2NR, where 
N is the number of columns in such a cluster. In order to 
find the entry energy due to such a cluster we will follow 
the Refs. [) 10, where the current and field distributions 
near the edge of the wedge-like surface crack where cal- 
culated. Consider thus a cut normal to the surface made 
of the chain of N columns. According to Refs. p|,[To| the 



current near the tip of a thin crack grows as jy^2NR/x 
when distance x from the tip of the crack (the last col- 
umn) falls into interval R < x <C 2NR, and saturates 
at jy2N for x < R. Here j is the Meissner current far 
away from the crack. The suppression of the barrier de- 
pends on the relation between the position of the energy 
maximum xq and R. 
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FIG. 9. The gate at the surface created by an improbable 
event when TV columns form a chain. Such gates dominate 
flux penetration at low temperatures. 

Consider the case x > R first. Making use of the 
results of Refs. ^,^0[ we write the energy of vortex at 
distance x from the tip of the crack ast3 



E(x) = ^jVWRx + se In 7. 

Accordingly, the barrier for this gate configuration (the 
maximum value of E(x)) is 



U = Uq — seq In 



QpjNR 
ce 



(10) 



where the bare barrier Uq in real units is given Uq — 
ssohi(jdp/j) and jd P is the depairing current. The en- 
ergy E(x) achieves its maximum value at the point 
xq = 2 (c£o/$oj) 2 / (NR), which has to satisfy conditions 
R < xq < NR. 

In order to estimate the effective penetration rate one 
has to determine the size of the clusters giving maximal 
contribution to vortex entry. The probability to find a 
chain of TV columns in a row scales as (TrR 2 rid) N with rid 
being the column concentration. Therefore the contri- 
bution from such chains to the penetration rate is given 
by 



H ex P 

N 



-TV In 



ceo 



1 £fo ln $ojN R\ 
TrR 2 rid T ceo 



)■ 



(u) 



and becomes maximal at the optimal TV 



TV = 



se 



> max 



Tln{l/irR 2 n d ) 
Now the condition xq > R reduces to 



2ce 



3>ojR 



T > 



se 



®ojR 

ln(l/irR 2 n d ) \ 2ce 



In this "high-temperature" regime the change of the sur- 
face barrier is: 



SU w — seo In 

= — seo m 



s<P jR 



2cTln(l/7ri? 2 n d ) 

sR<S> H 
8n\T \n(l /nR 2 n d ) 



where, again, j = cH/AttX. 

The effective penetration barrier is determined by the 
relation: 



with Ht 



H C H T n 
se a In 2 = CT 

s7?$n which gives 



8ttAT ln(l/rr R?n d ) ' 



H p = \J H C H T exp 



2CT\ 
se J 



(12) 



(13) 



Consider now the regime xq < R, corresponding to 
low temperatures. In this case the energy of the vortex 
is given by 

E(x) = -—jV2Nx + se In %. 

c 4 

The position of the energy maximum and the barrier are 
given by 



Xq 



ceo 



$oiV2/V 



U = seo m —? 



„C/ o _f££ ln7V 
$oOV2TV 2 



Following the same route one arrives at the optimal clus- 
ter length as 



TV 



se 



> 1, 



2Tln(l/7ri? 2 n d ) 
and, accordingly, the change of the barrier is 



SU 



se 



hi 



se 



2 2Tln(l/7ri? 2 n £ ;) 



Contrary to the case of individual columns, the cluster 
gate exhibits the temperature dependence: the depres- 
sion of the barrier decreases logarithmically with grow- 
ing temperature. The change of the effective penetration 
field is determined by 



, H c se 

seo m T7 77~ m 



se 



H p 2 2Tln(l/7ri? 2 n d ) 



CT 



or 



H„ 



l2T\n(l/7rR 2 n d ) 

se 



H, 



i>o 



(14) 



This formula gives the effective penetration fields at low 
temperatures. It is probably more relevant to experimen- 
tal situation than the "high-temperature" result jlgj). 
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We see that long clusters of columnar defects forming 
the crack-like configurations further (as compared to the 
effect of an individual column) suppress the surface bar- 
rier and serve as a very effective vortex gates into the 
sample. The collective mechanism wins over the single- 
column mechanism at low temperatures. Comparing Eq. 
( |l4| ) with the single-column result we conclude that the 
at high magnetic fields the crossover between the single 
column and collective regimes is expected at the crossover 
temperature T cr 

Tcr ~ 81n(l/7r°R2n d ) (15) 

For the typical parameters of the compound 
Bi 2 Sr 2 CaCu20 l: at rid = 5 ' 10 10 cm~ 2 (correspond- 
ing to the matching field of 1 tesla) and R = 35A the 
crossover is expected at w 20 K. 

Our calculations suggest that in the regime of the 
thermally-activated pancake penetration through the 
surface barrier, the penetration field is roughly two times 
smaller than the penetration field of unirradiated samples 
at temperatures above the crossover temperature ( |15| ) 
and decreases according to Eq. at lower tempera- 
tures. This is consistent with the recent experiment.! 7 ] 

We have investigated the influence of columnar defects 
in layered superconductors on the thermally activated 
penetration of pancake vortices through the surface bar- 
rier. Depending on the position of an isolated column 
the effective barrier is determined either by the vortex 
hopping from the surface to the column or by the de- 
tachment of the vortex from the column to the bulk. For 
a given external field there exists an optimum location 
of the column for which the barriers the both stages are 
equal and the reduction of the effective penetration bar- 
rier is maximal. Formation of long clusters of columnar 
defects thrusting the surface offers the most convenient 
gates for the vortex entry, the effect of cluster-induced 
depression of the surface barrier decreasing with tem- 
perature. Penetration through the clusters always dom- 
inates at low temperatures. It would be very interesting 
to investigate experimentally temperature dependencies 
of the surface barriers and penetration rates in order to 
reveal the role of clusters. The proposed mechanism of 
the disorder-assisted surface creep is very general and can 
be extended to the point disorder containing samples. 

We thank C. J. van der Beek, M. Konczykowski, and 
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ported by the U.S. DOE, Office of Science, under contract 
# W-31-I09-ENG-38. 



1 C. J. van der Beek, M. Konczykowski, R. J. Drost, P. H. 
Kes, N. Chikumoto, and, S. BoufTard, Phys. Rev. B 61, 
4259 (2000). 

2 V. N. Kopylov, A. E. Koshelev, I. F. Schegolev, and T. G. 
Togonidze, Physica C 170, 291 (1990). 

3 M. Konczykowski, L. Burlachkov, Y. Yeshurun, and F. 
Holtzberg, Phys. Rev. B 43, 13707 (1991). 

4 R. G. Mints and I. B. Snapiro, Phys. Rev. B, 47, 3273 
(1993). 

5 E. Zeldov, A. I. Larkin, V. B. Geshkenbein, M. Kon- 
czykowski, D. Majer, B. Khaykovich, V. M. Vinokur, and 

H. Shtrikman, Phys. Rev. Lett. 73, 10 (1994). 

6 L. Burlachkov, V. B. Greshkenbein, A. E. Koshelev, A. I. 
Larkin, V. M. Vinokur, Phys. Rev. B 50, 16770 (1994). 

7 J. K. Gregory, M. S. James, S. J. Bending, C. J. van der 
Beek, and M. Konczykowski, to be published 

8 I. M. Lifshitz and V. Ya. Kirpichenkov, Zh. Eksp. Teor. 
Fiz. 77, 989 (1979) (Sov. Phys. JETP 50, 499 (1979)) 

9 A. I. Buzdin and M. Daumens, Physica C 294, 257 (1998). 

10 A. Yu. Aladyshkin, A. S. Mel'nikov, I. A. Shereshevsky, and 

I. D. Tokman, Preprint, 2001, to be published in Physica 
C. 

11 D. Yu. Vololasov, Phys. Rev. B, 62, 8691 (2000). 

12 In layered superconductors with very weak interlayer 
Josephson coupling the image technique gives the energy of 
the point vortex near the surface with accuracy s/A. The 
exact solution and justification of the image approach has 
been given in Ref. A. Buzdin and D. Feinberg, Phys. Lett. 
A, 167, 89 (1992). Exact field distribution for an isolated 
pancake vortex has been calculated in Ref. R. G. Mints, 
I. B. Snapiro and E. H. Brandt, Phys. Rev. B, 54, 9458 
(1996). 

13 G. S. Mkrtchyan and V. V. Shmidt, Sov. Phys. JETP 34, 
195 (1972); A. Buzdin and D. Feinberg, Physica C 235- 
240, 2755 (1994), Physica C 236, 303 (1996). 

14 Electrostatic analogy for the problem of vortex interaction 
with defects has been extensively explored recently in the 
paper A. Buzdin and M. Daumens, Physica C, 332, 108 
(2000). 

15 Note that the column radius R appears in this equation 
only because our unit of current contains R. In real units 
this barrier, of course, does not depend on R. 

16 Because we only have qualitative estimates in this section, 
we use physical units from the very beginning. 



8 



